by GARTH GADDY 



How to Write a 
PID Algorithm 

The tricky part is trying to build 
a PID that operates smoothly in 
all modes of operation. 



When you hear of 
people you know 
having an 
angiogram or stents 
put in a collapsed 
artery, you realize these procedures 
can involve a whole variety of medical 
catheters of differing sizes. These 
catheters are manufactured by the 
process of polymer extrusion. Quality 
control of the catheters demand a pre- 
cise outside diameter with tolerances 
of +/-0.0005in. One method of main- 
taining an accurate outside diameter is 
to maintain accurate polymer pressure 
inside the extruder. The need for pre- 
cise pressure control is what brought 
me to the task of writing my own PID 
(proportional-integral-differential) 
algorithm. Those of you who have 
worked with PID controllers know that 
a wide variety of algorithms exist in 
the industry. The intent of this article is 
not to cover all PID algorithms but to 
show one implementation of a "home- 
made" PID and how I managed some 
of the trickier aspects of PIDs, such as 
bumpless transfer and setpoint 
changes. 

The extrusion systems I work on use 
a Delta Tau Data Systems PMAC (pro- 
grammable multi-axis controller) to 
control a servo motor that drives the 
extrusion screw that in turn provides 
the polymer pressure. In order to con- 
trol pressure, the PMAC AID convert- 
er reads the melt pressure and then 
passes data to my pressure PID whose 
output is then cascaded into the set- 
point of the extruder screw servo 
motor control loop. My PID is written 
in the PMAC motion control language, 
which is similar to Basic, and the code 



examples can be easily converted to C, 
C++, or other languages. 

Up until writing this algorithm, I 
had always depended on PIDs supplied 
with the control hardware/software 
package I was using. Whether I was 
working on a PLC, a motion controller, 
or a PC-based control system, I was 
only required to insert the PID into the 
control logic, run the program, and 
tune the loop. As I said earlier, this is 
an extruder pressure control loop 
whose output is cascaded into the set- 
point of a screw servo loop. One option 
I had was to employ an unused axis of 
the 8-axis motion controller for this 
pressure loop and cascade its output 
into the screw loop. This option would 
have wasted the motion control axis 
which might be needed later, so I 
decided to write my own loop and 
embed it in the screw servo motion 
program. Each time the servo motion 
program executed, the pressure loop 
would also execute and deliver a new 
screw loop setpoint. 

I started researching the subject and 
found many books on the subject of 
PIDs. Most were heavy on theory and 
either light or weightless on practical 
"show-me-the-code" kind of informa- 
tion I was looking for. Control system 
integrators I talked to knew how to 
tune one but not how to write one. It 
was black magic. I then found a book 
on the basics called Repl-Time 
Computer Control: An Introduction 
(Stuart Bennett, Prentice Hall 
International, Hertfordshire, UK, 
1994). The book covered writing a 
generic PID loop and also included 
some methods of bumpless transfer. I 
put together a basic PID and had the 
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Until writing this 
algorithm, I had 
always depended 
on PIDs supplied 
with the control 
hardware/ 
software package 
I was using. 



folks at Delta Tau check my algorithm. 
They gave it a thumbs up but men- 
tioned that I might have difficulty 
going from manual to auto and back 
without further code. I was so excited 
to have something that was close to 
implementation that I said I would fig- 
ire it out and headed for the extrusion 
lab. As soon as I ran some test numbers 
through the algorithm I knew that they 
weren't kidding. More code was need- 
ed. I didn't like the bumpless transfer 
methods in the book, and after banging 
my head against the extruder for a few 
hours, I came up with the method 
described in the examples. 

Many differing PID algorithms are 
used in industry. Each PLC manufac- 
turer that offers a PID adds a custom 
twist to the algorithm. Each PC-based 
control system that has a PID also adds 
a custom twist. The only exception is 
the classic ISA PID algorithm that 
conforms to a standard. For the pur- 
poses of my examples below, I will be 
referring to the independent or parallel 
PID algorithm (which I used). See the 
sidebar "Abbreviations" for a few def- 
initions of the abbreviations I use in the 
following text and code fragments. 

ERROR 

As its name implies, a PID con- 
sists of three components or 
terms: proportional, integral, 
and derivative. The proportional term 
is a product of error and proportional 
gain. Error is calculated by finding the 



difference between the setpoint and the 
process variable. This can be done by 
either: 

Error = SP - PV 
or: 

Error = PV - SP 

The method you select will deter- 
mine whether your control loop is a for- 
ward or reverse acting PID (some call 
this either inc/inc or inc/dec). In other 
words, as the output of the loop increas- 
es, does the process variable go up or 
go down? How your process responds 
will determine which method you 
choose. For example, I used the first 
method above of calculating error in 
my polymer pressure loop PID. An 
increase in the output of that loop 
(servo velocity command) turns the 
screw faster, which in turn yields a 
higher pressure. Another way to look at 
it is if there is positive error in Equation 
2, the pressure is lower than the set- 
point. Positive error will always yield 

ABBREVIATIONS 

■ Setpoint (SP) 

■ Process Variable(PV) 

■ Proportional Gain (KP) 

■ Integral Gain (Kl) 

■ Derivative Gain (KD) 

■ PV High Limit - PV Low Limit 
(PV_Span) 

■ CV High Limit - CV Low Limit 
(CV_Span) 

■ PID or Manual Flag (Auto_Mode) 

■ Error (Error) 

■ Sum of the Error (Sum_Error) 

■ Change in Process Variable Since Last 
PID Execution (d_PV) 

■ Change in Error Since Last PID 
Execution (d_Error) 

■ Proportional Term (PTerm) 

■ Integral Term (ITerm) 

■ Derivative Term (DTerm) 

■ Process Variable on Last PID Execution 
(Last_PV) 

■ Error on Last PID Execution (Last_Error) 

■ Error when manual to auto 
(lnitial_Error) 

■ Last Setpoint when setpoint changes 
(Last_SP) 

■ Controlled Variable (CV) 

■ Manual RPM Setpoint 
(Manual_Command) 

■ PID Execution Stage 
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an increase in the loop output (we'll see 
that later), which, in this case, will turn 
the screw faster and increase the pres- 
sure, thus reducing the error. If I had 
chosen the second method and there 
had been any error, the loop would 
have sent the screw velocity in the 
wrong direction (faster for high pres- 
sure or slower for low pressure) and the 
pressure would have "run away" in the 
direction of the error. 

An example of using the method in 
Equation 2 is a tank level control using 
a proportional valve to let liquid out as 
an amount of liquid is being pumped 
in. In this case, the PV is the actual 
level and the SP is the desired level. 
The output of the loop controls the per- 
cent the discharge valve is open, from 
0% (fully closed) to 100% (fully open). 
After a bit of thought it is apparent that 
Equation 2 is the best choice. Let's do 
a what-if to make sure. If the tank were 
too low, it would cause negative error, 



FIGURE 1 

Error over time. 



error 




time 




which would then reduce the controller 
output, closing the valve, and bring the 
level back up. If Equation 1 were cho- 
sen, a low tank would cause positive 
error, open the valve more, and empty 
the tank. 

PROPORTIONAL TERM 

The proportional term is so 
named because its output is 
proportional to error. It 
responds to instantaneous error, error 

\ 



that is calculated on that one execution 
of the PID. If you had a P-only con- 
troller, with a Kp greater than zero, the 
controller's output will increase as the 
error increases and decrease as error 
decreases. The typical (simplified) 
code fragment for a P-only controller 



is: 



Pterm = Kp * Error 
CV = Pterm 
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As you can also see, the higher the 
Kp the higher the output for a given 
amount of error. That's why a loop 
with a high proportional gain is more 
responsive to error than one with a 
respectively lower gain. 

Those of you with a perceptive eye 
might have noticed something that did- 
n't seem like it would work: units. 
How do you handle dissimilar input 
and output units? My pressure PID SP 
and PV units are in psi, while the CV is 
in RPM. The process of converting the 
units is called normalization and is 
handled by setting up a ratio of the 
input and output spans, and multiply- 
ing the result by the variable you want 
to convert (error in this case). The final 
Pterm code fragment is as follows: 

Error = (SP - PV) * (CV.Span / PV.Span) 



or: 



Error = (PV - SP) * (CV.Span / PV.Span) 



;Then 

Pterm = KP * Error 
CV = Pterm 

INTEGRAL TERM 

The integral term is so named 
because its output is the integral 
of the error with respect to time. 
Its job is to remove long term or steady 
state error that is not removed by the 
Proportional Term. As you may 
remember from calculus, the first inte- 
gral of an equation of a line is equal to 
the area under the line. In Figure 1, the 
line shows error over time and the area 
under the line is the integral of the 
error. We calculate the integral by 
adding the error to itself each time the 
PID executes. Integral term code frag- 
ment in an integral-only controller 
looks like this: 

Error = (SP - PV) * (CV.Span / PV.Span) 
or: 



Error = (PV - SP) * (CV.Span / PV.Span) 

;Then 

Sutn.Error = Sum.Error + Error 
Item = KI * Sum.Error 
CV = Item 

If an error persisted that was not cor- 
rected by the proportional term and the 
KI is greater than zero, the integral 
term would begin to grow larger and 
larger eventually correcting or "reset- 
ting" the long term error. 

DERIVATIVE TERM 

The derivative term is meant to 
slow down changes in the 
process being controlled. This 
slowing is done by calculating the 
changes in the process or the "deriva- 
tive" of the process with respect to 
time. Two indicators of process change 
are available, PV and error. Derivative 
Term code fragment is as follows: 

;PV Method 
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LISTING 1 

Screw servo motion program. 



01 

02 
03 
04 
OS 
06 
07 
08 
09 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 



While (True) 

If (Huto.Hode = 1 ) 

If (First.PID.Execution = OR First.PID.Execution = 1) 
Error = (SP - PV) * (CV.Span / PV.Span) 
Sumjrror = Sum.Error ♦ Error 
If (First.PID.Execution = 0) 

Sum.Error = HanuaLCommand / KI 
First.PID.Execution = 1 
Initialjrror = Error 



Endif 
Pterm ■ 
Item ; 
Dterm ; 



Else 



KI * Sum.Error 

: 

If (Initial.Error > and Error < 0) 

First.PID.Execution * 2 
Last.PV = PV 

Endif 

If (InitiaLError < and Error > 0) 

First.PID.Execution = 2 
Last.PV = PV 

Endif 

Last.SP = SP 

If (SP = Last.SP) 

Error = (SP - PV) * (CV.Span / PV.Span) 

Pterm = KP * Error 

Sum.Error = Sum.Error + Error 

Iterm = KI * Sum.Error 

d.PV = (Last.PV - PV) * (CV.Span / PV.Span) 

Last.PV = PV 

Dterm = KD * d.PV 

Else 

Error = (SP - PV) * (CV.Span / PV.Span) 
Sum.Error ■ Sum.Error + Error 
Pterm = 

Iterm = KI * Sum.Error 
Dterm = 

If (SP > Last.SP and PV > SP) 

Last.SP = SP 
Last.PV = PV 

Endif 

If (SP < Last.SP and PV < SP) 

Last.SP = SP 
Last.PV = PV 

Endif 

Endif 




Endif 

CV ■ Pterm + Iterm + Dterm 



Else 



CV ■ Manual_Command 
First.PID.Execution = 



Endif 
EndWhile 



Error = (SP - PV) * (CV.Span/PV.Span) 
d.PV = (Last.PV - PV) * (CV.Span/PV.Span) 
Dterm = KD * d.PV 
Last.PV = PV 
CV = Dterm 

; Error Method 

Error = (SP - PV) * (CV.Span/PV.Span) 
d.Error = Last.Error - Error 
Dterm = KD * d.Error 
Last.Error = Error 
CV = Dterm 

As you can see, if the PV method 
used, the normalization code must 
added in. If the error method is use 
normalization has already been do 
in the error calculation. The output 
the derivative term will be large if K 
is greater than zero and there is a lar 
change in either PV or error. Changii 
the derivative gain will change ti 
effect that the Dterm has on the fin 
output. 

PUTTING IT ALL TOGETHER 

In a complete PID algorithm, ; 
three terms are put together. Ti 
code fragment that follows is a fc 
ward acting, PV based Dterm Pl6: 

Error = (SP - PV) * (CV.Span / PV.Span) 

Pterm = KP * Error 

Sum.Error = Sum.Error + Error 

Iterm = KI * Sum.Error 

d.PV = (Last.PV - PV) * (CV.Span / PV.Spai 

Last.PV = PV 

Dterm = KD * d.PV 

CV = Pterm + Iterm + Dterm 

So that's it. How to write a PID. ( 
to it and good luck. But wait, will th 
thing really work? You know, in tl 
real world? Well, yes and no. This PI 
would work but not transition from 
manually controlled process to a PII 
controlled process very gracefully, 
would also not handle setpoint chang 
very well. So let's take a look at tho 
two issues. Because the remainii 
code will have more features, refer 
Listing 1 in order to follow along in t 
text. 
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BUMPLESS TRANSFER AND SETPOINT 
CHANGES 

Bumpless transfer is the act of 
taking a process from manual- 
ly-controlled to PID-controlled 
without any noticeable upset or 
"bump" in the process. I will use my 
extruder pressure control PID for the 
next examples in order to drive home 
this concept. 

Let's say my extruder screw is run- 
ning manually controlled at 15 RPM. If 
I were to use the above code alone and 
put the screw in auto (PID loop on), the 
output of the loop on the first execution 
would go from 15 RPM in manual to 
something barely above zero, based on 
the pressure error at the instant the loop 
executed. The output would then head 
off towards the setpoint at a rate depen- 
dent upon gain values. Not very bump- 
less. In order to make the transition 
bumpless, I need to somehow have the 
output of the PID loop equal the man- 
ual output at the instant of manual to 
auto switchover. I do this by "pre-load- 
ing" the Iterm with the manual output 
command on the first execution of the 
PID. The loop is then kept in an inte- 
gral-only mode until the PV meets the 
SP, and then it is switched into PID 
mode. 

Smoothing out setpoint changes will 
allow the operator to make large set- 
point changes without upsetting the 
process. This outcome is often 
achieved with a "ramped setpoint." I 
chose a method which resulted in less 
code by switching the controller to an 
integral-only mode until the process 
reaches the setpoint. The controller is 
then switched back into PID mode. 

Let me say a few words about 
Listing 1: 

■ The PID is set in an infinite loop 
between Line 01 and Line 53 

■ It is then broken down into two 
modes, auto and manual, which are 
determined by the If-Then-Else 
statement starting on Line 02, Else 
at Line 49, and ending at Line 52 

■ The Auto section is from Line 03 to 
Line 48 and has two subsections 



■ The Auto subsection 1 is the bump- 
less transfer section from Line 04 to 
Line 22 

■ The Auto subsection 2 is the normal 
operation section from Line 24 to 
Line 46 

■ Within the normal operation sec- 
tion, there is one more set of sub- 
sections 

■ Normal operation subsection 1 is 
PID control from Line 25 to Line 3 1 

■ Normal operation subsection 2 is 
for setpoint changes from Line 33 to 
Line 45 

In order to understand how the 
Listing 1 PID works, let's walk 
through a typical scenario of startup 
and transition into Auto or PID mode. 
Prior to start, the operator makes sure 
the machine is set in manual mode and 
then enables the screw servo loop. This 
action executes the Listing 1 code 
within the infinite loop (Lines 01 to 
53). The operator then enters a 
Manual.Command (RPM) on the operator 
interface and the screw begins turning 
at that RPM. The only statements that 
are performing any servo control at 
this point are Lines 50 and 51. These, 
lines keep the servo at its manual set- 
point and reset the First.PID.Execution 
flag. When polymer pressure is up to 
desired level, the operator places the 
screw in Auto (PID mode). From here 
on the only control code executed will 
be between Lines 02 and 49. On first 
scan in PID mode, Line 07 sets up the 
Iterm "pre-load", Line 08 sets the 
First.PID.Execution flag so these three 
lines of code will only run once, and 
Line 09 sets an Initial_Error variable 
to the error. This will be used in con- 
junction with Lines 14 to 22 to gently 
bring the extruder pressure up to set- 
point. The next three lines of code set 
the three PID terms and you'll, notice 
that the Pterm and Dterm are set to 
zero. This leaves an I-only controller 
whose Iterm will slowly (gently, 
bumplessly) get larger or smaller 
depending on the error until it forces 
the output to a point where the PV 
passes the SP. This event is picked up 
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by either Lines 14 to 17 or 18 to 2 
depending on the sign of th 
Initial_Error and the controller is the 
shifted into PID mode from Lines 25 1 
31. Before we leave the bumples 
transfer code explanation, note Line 
16 and 20 — Last.PV is set equal to P\ 
This is necessary to assure that on fir; 
calculation of the Dterm, its output i 
stable. Otherwise a severe "bump 
may result on transition. 

Lines 25 to 31 should be sell 
explanatory by now for they are th 
standard, no-frills PID. If the operato 
should select a new pressure setpoir 
(SP) code from Lines 33 to 45 execute 
until the PV passes the new SP. You'i 
notice that this code looks similar t. 
the bumpless code mentioned' above 
The code uses the same technique o 
turning the PID into a straight I-con 
troller, which is more gentle abov 
seeking a large setpoint change than ; 
full PID. This section of code is exe 
cuted until either Line 38 or 4. 
(depending on the relationshi] 
between SP and Last.SP) see that th 
SP has been passed by the PV. It i 
then shifted back to standard PID. 

SMOOTH OPERATOR 

As you can see, PID code alon 
is fairly simple. The trickie 
part comes in when trying t< 
build a PID that operates smoothly ii 
all modes of operation. Many way 
exist to accomplish PID control, bump 
less transfer, and handle setpoin 
changes. Other methods that could am 
probably should be tried are variabl 
gains and ramped setpoint. Eac) 
method must be weighed against th' 
type of PID algorithm used. Wha 
worked on a parallel algorithm ma; 
not work well on an ISA PID. WSQ 
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